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We present molecular simulation data for viscosity, self-diffusivity, and the local structural ordering of (i) a 
hard-sphere fluid and (ii) a square-well fluid with short-range attractions. The latter fluid exhibits a region of 
dynamic anomalies in its phase diagram, where its mobility increases upon isochoric cooling, which is found 
to be a subset of a larger region of structural anomalies, in which its pair correlations strengthen upon isochoric 
heating. This "cascade of anomalies" qualitatively resembles that found in recent simulations of liquid water. 
The results for the hard-sphere and square-well systems also show that the breakdown of the Stokes-Einstein 
relation upon supercooling occurs for conditions where viscosity and self-diffusivity develop different couplings 
to the degree of pairwise structural ordering of the liquid. We discuss how these couplings reflect dynamic 
heterogeneities. Finally, we note that the simulation data suggests how repulsive and attractive glasses may 
generally be characterized by two distinct levels of short-range structural order. 



I. INTRODUCTION 

Colloidal fluids play a central role both in technology and 
in the study of condensed matter. Regarding the latter, sus- 
pensions of colloids are interesting model systems because 
they can behave collectively in ways that are similar to atomic 
and molecular liquids while simultaneously being both large 
and slow enough to allow experimental measurements of their 
real-space structure and dynamics. Despite the similarities be- 
tween the behaviors of colloids, atoms, and molecules, there 
also are important, qualitative differences^^ For example, 
the effective interparticle attractions between colloids can be 
"tuned" to be short-ranged relative to the particle diameter cr, 
unlike the dispersion attractions in atomic and molecular sys- 
tems that decay more gradually with interparticle separation 
r, i.e., as (cr/r)^^. 

Short-ranged attractive (SRA) interactions have nontrivial 
implications for the static structural, equilibrium thermody- 
namic, and dynamic behaviors of colloidal suspensions. For 
instance, whereas cooling a simple atomic liquid generally 
slows down its dynamic processes, several recent studies have 
demonstrated that reducing the temperature (or increasing 
the attractive interactions) of an SRA fluid can have a non- 
monotonic effect on its mobility. In fact, concentrated SRA 
fluids can vitrify not only upon cooling, forming an "attrac- 
tive" glass or gel, but also upon heating, forming a "repul- 
sive" or hard-sphere (HS) glass. Since the structural 
and mechanical properties of these glassy states can be quite 
different, there is a broad interest in understanding the prop- 
erties of the precursor supercooled fluids from which they are 
formed. 

In this article, we explore whether the unusual effect that 
temperature has on the dynamics of SRA fluids reflects a more 
general connection between the static structure and the dy- 
namics of condensed phases. In doing so, we find it instruc- 



tive to view the behavior of SRA fluids in the context of an- 
other well-studied system with anomalous dynamical trends, 
liquid water. Cold water behaves differently than simple flu- 
ids in that its mobility increases upon isothermal compres- 
sion over a broad range of conditions. Interestingly, the state 
points where water shows this behavior is a subset of a larger 
region on its phase diagram where its local structural order 
anomalously decreases upon compression. In other words, 
there is a "cascade of anomalies"'' where the pressure-induced 
disordering of liquid water emerges at a lower density, and 
ultimately disappears at a higher density, when compared to 
the pressure-induced increase in its mobility. It has been re- 
cently arguediiii^iil that this behavior follows from the fact 
that water approximately obeys a scaling relationship between 
its self-diffusivity and the structural order parameter —52/^6, 
where S2 is the contribution to the fluid's excess entropy due 
to its static oxygen-oxygen pair correlations, and fee is the 
Boltzmann constant. 

Here, we use molecular dynamics simulations to systemati- 
cally investigate the relationships between static structural or- 
der (—82/ kB), self-diffusivity (D), and viscosity (77) for both 
a HS fluid and an SRA fluid. One of our main aims is to 
understand whether the latter shows a cascade of anomalies 
similar to that of liquid water. Specifically, we are interested 
in whether the region on the phase diagram of an SRA fluid 
where it becomes more mobile (higher D, lower rj) upon cool- 
ing is a subset of a larger set of conditions where the fluid 
becomes more ordered (higher —52/^8) upon heating. A sec- 
ond, and related, goal of this study is to explore whether the 
relationships between D, 77, and —82/ can provide generic 
insights into the breakdown of the Stokes-Einstein (SE) rela- 
tion (Dri/T « constant) in deeply supercooled liquids and 
also into the structural properties of repulsive and attractive 
glasses. 

The organization of the paper is as follows. In Section HIl 
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we describe the two model systems examined in this study 
and also the simulation methods used to carry out the investi- 
gation. Then, in SectionlnD we present our simulation results 
and discuss their relevance for understanding the connection 
between structural order and mobility in HS and SRA fluids. 
Finally, in Section lTVl we present some concluding remarks. 



11. MODELING AND SIMULATION 

We focus on two model systems: a fluid of HS particles and 
a fluid of square-well particles with short-range attractions, 
the latter of which we denote the SW-SRA fluid. The inter- 
particle potential V(ri2) between particles 1 and 2 for these 
systems is given by 



using the expression 

- S2 = 1 1 {3(r)ln5(r) - [^(r) - l]}dr, (2) 

where p = N/V is the number density, and g (r) is the average 
pair correlation function. Note that —S2 — for an ideal gas, 
and — S2 oo for a perfect crystalline lattice. Thus, as has 
been discussed at length elsewhere one can view — S2 
as a scalar measure of the amount of pair-wise translational 
order of the system. 



III. RESULTS AND DISCUSSION 



A. Transport and Structural Properties 



oo ri2 < cri2, 
V(ri2) = < cri2 < ri2 < Acri2, 
ri2 > Xai2- 



(1) 



where ri2 is the distance between the particle centers, cti2 = 
(cti + o"2)/2 is the interparticle diameter, and the parame- 
ters e and A set the magnitude and range of the interparticle 
attraction, respectively. For the HS fluid, one has A = 1, 
and thus there are no interparticle attractions. We choose 
A = 1.03 for the SW-SRA fluid, a range similar to that of 
other SW-SRA fluids known to exhibit anomalous dynami- 
cal behavior** '**. In order to avoid crystallization in our study, 
and thus allow study of the supercooled fluid states, we have 
drawn the individual particle diameters di of both systems 
from a Gaussian distribution with an average of cr and stan- 
dard deviation of s = 0.1a. For practical reasons, we trun- 
cated this distribution so that all particle diameters lie in the 
range cr — 3s < cTj < a + 3s. We have implicitly non- 
dimensionalized all reported quantities in this investigation 
by appropriate combinations of the characteristic length scale, 
lc = cr and a model-dependent time scale t^. For the HS fluid, 
the characteristic time scale Tc is cast in terms of the temper- 
ature Tc ^ma^ jk^T, where m is the particle mass, while 
for the SW-SRA fluid it is defi ned in te rms of the attractive 
strength of the potential, Tc = ^Jma^ je. 

To explore the behavior of these model systems, molec- 
ular dynamics simulations were performed using an event- 
driven algorithm'^ in the microcanonical ensemble. For all 
runs, iV = 1000 particles were simulated in a cubic simu- 
lation cell of volume V with periodic boundary conditions. 
Particle packing fractions 0c = "^^^Z^t I ^^^^ range 
0.35 < 0c < 0.6 were investigated for both model systems, 
and temperatures T in the range 0.3 < T < 2.0 were exam- 
ined for the SW-SRA fluid. For each thermodynamic state 
point, between three and ten independent simulations were 
performed in order to estimate the errors in the transport co- 
efficients. Self-diffusivities D were calculated by fitting the 
long time (t 1) mean-squared displacements to the Einstein 
relation for self-diffusivity ((5r^) = 6Dt. Zero-shear viscosi- 
ties were calculated using the impulse limit of the Einstein 
formula.'^ The structural order parameter — S2 was computed 



We now examine some of the state-dependent transport and 
structural properties of the HS and SW-SRA fluids. For the 
HS fluid, the quantities D, rj, and — S2 (non-dimensionalized 
as outlined in Section|ll]i depend only on the packing fraction 
<j>c, whereas the corresponding dimensionless properties for 
the SW-SRA fluid depend on both packing fraction (pc and 
temperature T. 
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FIG. 1 : Transport properties of the HS fluid described in the text as a 
function of packing fraction (pc. (a) self-diffusivity D, (b) viscosity 
rj, and (c) SE relationship Dij/T. The horizontal dashed line in (c) 
indicates {2it)^^, the expected value of the SE relation in the slip 
limit, and the vertical line denotes the point at which which Dij — 
1.2/(27r). In this work, we use this simple heuristic to identify the 
breakdown of the SE relation. 
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FIG. 2: Transport properties of the SW-SRA fluid described 
in the text as a function of packing fraction <^c and reciprocal 
temperatureT^^: (a) self-diffusivity D, (b) viscosity tq, and (c) 
SE relationship D-q/T. The horizontal dashed line in (c) indicates 
(27r)^^, the expected value of the SE relation in the slip limit. 



The transport properties of the HS fluid are displayed in 
Fig.[T] As expected, with increasing packing fraction (pc, the 
self-diffusivity D monotonically decreases, while the viscos- 
ity rj monotonically increases. The changes in both of these 
properties become more pronounced for 0c > 0.55, condi- 
tions for which D{(j>f.) can be accurately fitted using either 
Vogel-Fulcher or power-law functional forms'^ with a pre- 
dicted divergence near (pc ^ 0.6. 

Figure [3c) illustrates that the SE relation in the slip limit, 
Drj w {2tt)~^, is approximately obeyed by this HS fluid for 
(f)c < 0.55. However, for higher packing fractions, large pos- 
itive deviations from the slip limit of the SE relation become 
apparent."^ This type of SE "breakdown" has been observed 
in studies of several glass-forming fluids in their deeply su- 
percooled states. '''•^'^'^''^^ 

We note that, as observed by previous investigators^i^ the 
SE relation is not strictly obeyed by the HS fluid even at lower 
packing fractions. For example, the product Drj varies by 
approximately 30% as the packing fraction is increased from 
(j)c = 0.35 to 0.55. However, the breakdown of the SE rela- 
tion is typically understood to occur when the product begins 
to exhibit a pronounced positive deviation from the slip value. 
In this work, we identify the breakdown of the SE relation by 
the heuristic criterion, Drj > 1.2/(27r). As can be seen from 
Figure[Tfc), this threshold is crossed at ipc ~ 0.55 for the HS 
fluid investigated here. 



The simulated transport properties of the SW-SRA fluid are 
displayed in Fig. |2] as a function of reciprocal temperature 
T^^ along isochores. The most striking feature of this plot 
is that D exhibits a maxima with inverse temperature, a be- 
havior that has also been observed in both experiments and 
computer simulations of other SRA fluids i^i'^i^'^'^'^i^^'^'*'^^ As 
discussed in Section|I] this trend becomes pronounced at high 
(pc, where the system can ultimately form either a repulsive 
glass by isochoric heating or an attractive glass or gel via iso- 
choric cooling. 

The behavior of the zero-shear viscosity 77 of the SW-SRA 
fluid, displayed in Fig. HJb), qualitatively mirrors that of its 
self-diffusivity. For a more quantitative comparison, the SE 
relationship Drj/T is plotted in Figure|2jc). Note that the slip 
limit of the SE relation is again approximately obeyed for the 
SW-SRA fluid over a broad range of T^^ and (pc- However, 
for all values of 4>c investigated here, the SE relation breaks 
down at sufficiently low T, as the attractive glass transition is 
approached. At high T, on the other hand, only the highest 
(pc isochore studied showed a significant breakdown of the SE 
relation. This asymmetry between high and low T behaviors 
has also been observed in other model SRA fluids. ■^'ii^iSl It 
simply reflects the fact that one must reach a relatively high 
value of (pc in these systems to achieve the level of frustration 
required to form a repulsive glass, whereas interparticle at- 
tractions can induce formation of the attractive glass at much 
lower particle concentrations. 




FIG. 3: Structural properties of the HS fluid discussed in the text, (a) 
Translational structural order parameter —S2 versus packing fraction 
(j>c. (b) Radial distribution function g{r) for several packing frac- 
tions, (c) Cumulative order integral (see Eq. [3}. Arrows indicate 
increasing (j>c- 



4 



It is natural to wonder whether the decrease in mobihty that 
the HS and SW-SRA fluids experience upon approach to the 
glass transition is generally correlated with an increase in the 
amount of local structural order that they exhibit. To explore 
this issue, we first examine the behavior of the structural order 
parameter — S2 for the HS fluid. As can be seen in Fig. |3la), 
—S2 for this system monotonically increases with packing 
fraction, indicating a strengthening of the pair-wise interpar- 
ticle correlations. From the radial distribution functions g{r) 
shown in Fig.[3jb), it is also evident that these correlations cor- 
respond to the progressive development of well-defined coor- 
dination shells around the particles. The contributions of these 
shells to the translational order parameter — S2 become readily 
apparent when we investigate the following cumulative order 
integral, 

h,{r)^27rpf r'^g{r')lngir')-[g{r')-l]}dr'. (3) 

Note that one recovers —S2 from this integral in the large r 
limit, and thus Ig^ (r) quantifies the average amount of trans- 
lational ordering on length scales smaller than r surrounding a 
particle. In Fig.[3c), we observe nearly step-wise increases in 
{r) at the locations of the various coordination shells. It is 
also clear that increasing 0c of the HS fluid has two main ef- 
fects on its structural order. It strengthens the ordering within 
the individual coordination shells, and it uniformly increases 
the number of coordination shells (i.e., the range of order). 

The behavior of the structural order parameter — S2 for the 
SW-SRA fluid is displayed in Fig.|4ja) as a function of recip- 
rocal temperature T^^ along isochores. In accord with what 
might be expected based on the behavior of this fluid's trans- 
port coefficients, —S2{T) displays broad minima. At low T, 
heating the fluid decreases its structural order (similar to what 
happens in normal molecular fluids), but, at high T, heating 
anomalously increases its structural order 

To determine the origins of this trend, we focus on the be- 
haviors of the radial distribution function and the cumulative 
order integral for the SW-SRA fluid along the (pc = 0.55 iso- 
chore. As can be seen in Fig.|4jb)-(e), heating induces subtle 
changes in g{r) that result in nontrivial cumulative changes 
in the structural order In particular, heating the cold fluid 
(T < 0.4) results in a barely detectable decrease in the height 
of the first peak of g{r), and it has little effect on the other 
coordination shells. However, as is shown in Fig. |5c), even 
these subtle modifications to g{r) result in significant changes 
to the cumulative order integral of the fluid. In particular, 
heating the cold fluid decreases the total amount of pair-wise 
structural order, even though the range of the order remains 
essentially the same. This change is qualitatively consistent 
with the thermal weakening of transient clustered configura- 
tions of physically "bonded" particles observed in SRA fluids 
at lower temperature Si2i2^i2^ 

The effect of increasing temperature on the structure of the 
warm SW-SRA fluid (T > 0.4) is different [see Fig.Hd) and 
(e)]. In particular, heating continues to decrease the height of 
the first peak in g{r), but it also broadens the peak and en- 
hances the interparticle correlations associated with the other. 
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FIG. 4: Structural properties of the SW-SRA fluid, (a) Structural 
order parameter — S2 versus reciprocal temperature T^^. Symbols 
are the same as in Fig.|2l (b) Radial distribution function g{r) and 

(c) cumulative order integral /sal'") for <l)c ~ 0.55 and T < 0.4. 

(d) Radial distribution function g{r) and (e) cumulative order inte- 
gral (r) for 0c = 0.55 and T > 0.4. 



more distant, coordination shells. The net effect is an increase 
in both the total amount of translational structural order and 
its range. This structural change is due to the fact that heat- 
ing the warm fluid collapses the open channels of free vol- 
ume that form at intermediate T due to weak interparticle 
clustering, essentially jamming the particles into a less ef- 
ficient, and more correlated, packing arrangement. iii^fS 

It is clear from the results presented in this section that there 
is a qualitative (negative) correlation between pair-wise struc- 
tural order and mobility in both the HS and SW-SRA fluids. 
In the following two sections, we explore the extent to which 
this connection can be made quantitative. 



B. Connection between structure and mobility anomalies 

As discussed in Section UlI Al at high T and ipc, both the 
self-diffusivity and the structural order of the SW-SRA fluid 
behave in a manner that is anomalous when compared to sim- 
ple molecular fluids [see Figs. |2a) andUa)]. Whereas in- 
creasing T of a simple fluid generally increases its mobility, 
isochorically heating the warm SW-SRA fluid can result in 
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slower single-particle dynamics. This latter behavior is char- 
acterized by the following inequality: 



9D\ 



self-diffusivity anomaly. 



(4) 



The structural order of a simple molecular fluid, on the other 
hand, normally decreases when it is isochorically heated. 
Therefore, we denote conditions for which the following in- 
equality holds (i.e., order increases upon heating). 



d[- 



-S2\ 



dT 



> 0, structural anomaly, 



(5) 



as "structurally anomalous". 

The locations of the regions for self-diffusivity and struc- 
tural anomalies of the SW-SRA fluid in the T-ipc plane, as de- 
termined by numerical differentiation of the data in Figs.|2ta) 
andHfa), respectively, are displayed in Fig.|5ja). A schematic 
based on the data is provided in Fig. |5jb). The most strik- 
ing point is that the region of structural anomalies appears 
to completely envelop the region of self-diffusivity anoma- 
lies. In other words, the SW-SRA fluid exhibits unusual T- 
dependencies for its single-particle dynamics only for those 
state points where it also exibits unusual T-dependencies for 
its structural order. 

This "cascade" of structural and dynamic anomalies shown 
by the SW-SRA fluid is very similar to that observed origi- 
nally in simulations of the SPC/E model of water^ and later 
in simulations of other simpler models that also show water- 
like behavior. i ?-2B,29,30,3i Rgc^ij j-jjat cold water behaves dif- 
ferently from simple fluids over a wide range of conditions in 
that its mobility increases, while its structural order decreases, 
when it is isothermally compressed. Thus, the generic simi- 
larity between the behavior of water and the SW-SRA fluid 
shown in Fig.|5]is that, in both cases, the region on the phase 
diagram where mobility anomalies occur is a subset of the 
region where structural anomalies are found. It has been re- 
cently arguediiii^iil that this type of behavior for water fol- 
lows directly from the fact that liquid water approximately 
obeys a scaling relationship between its self-diffusivity and 
its translational structural order parameter — S2 over a wide 
range of temperature and density. In the next section, we test 
whether there exist similar scaling relations between —S2 and 
the transport coefficients of the HS and SW-SRA fluids. We 
also discuss how such relations might provide insights into the 
breakdown of the SB relation for these systems. 



Structure-property relations and the breakdown of 
Stokes-Einstein 
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FIG. 5: Conditions exhibiting self-diffusivity and structural anoma- 
lies, as well as breakdown of the SE relation [D-q/T > 1.2/(2%), 
see discussion in text] for the SW-SRA fluid in the T-(j)c plane, (a) 
Results from simulations. Large closed circles are state points where 
the SE relationship breaks down. Open circles represent the region 
of structural anomalies defined by Eq. lO. Small closed circles rep- 
resent the region of self-diffusivity anomalies defined by Eq. (b) 
Schematic representation of the data. The green shaded region (and 
area to its right) represents state points where the fluid is structurally 
anomalous. The blue shaded region represents state points exhibit- 
ing the self-diffusivity anomaly. Points to the right of the red curve 
show a breakdown of the SE relation. The gray region represents the 
repulsive and attractive glassy states. 



As alluded to above, recent molecular dynamics 

simulations", motivated by other earlier observations of model fluids in their equilibrium liquid states: 
Rosenfeld^^'-'^' and Dzugutov^'*, have demonstrated that the 

following simple relation is approximately obeyed by various D = Ad exp[i?£) S2], (6a) 
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where Ad and Bd are parameters which may depend on 
packing fraction (i.e., density), but not on T. Resuhs from 
other eariier theoretical studies, and considerations based on 
the Stokes-Einstein relation, suggest that a similar relation- 
ship should approximately hold for the zero-shear viscosity rj 
of these equilibrium fluids: 



1] ^ Ar, exp[B^ S2] 



(6b) 



where, again, A,j and S,, may depend on packing fraction 
(density) only. 

However, given that the SE relation breaks down as a Uq- 
uid is supercooled, it is apparent that Eq. ( l6al i and (l6bl i, with 
coefficients fit to higher temperature equilibrium fluid data, 
cannot also describe the transport coefficients in deeply super- 
cooled liquid states. Nonetheless, it has recently been shown 
that the functional form of Eq. (|6a] i can in fact approximately 
describe the isochoric self-diffusivity data of model super- 
cooled liquids over a broad range of temperatures if a differ- 
ent pair of parameters A'jj and B'j^ are adoptedjii In other 
words, there appears to be a crossover upon cooling where the 
self-diffusivities of fluids transition from being approximately 
described hy D = Ad cxp[i?/) S2] for equilibrium states to 
being approximately described hy D — A'j^ Gxp[_B^ §2] for 
supercooled conditions. Of course, even this latter relation 
must eventually fail for conditions very near the glass transi- 
tion where both D and 77^^ rapidly vanish, while S2 remains 
finiteiii We will return to this last point at the end of the sec- 
tion. 

The structure-property scalings discussed above suggest 
several interesting questions concerning the liquid state. For 
example, does the aforementioned crossover between scaling 
relations occur near the breakdown of the SE relation? Fur- 
thermore, does a similar crossover for the S2 dependence of 
r] occur upon supercooling? If so, does it coincide with the 
crossover point for Dl Below, we use our molecular dynam- 
ics simulation results to investigate these questions for the HS 
and SW-SRA fluids. The idea is that the new information that 
we gain about how D and 77 couple to pair structure should 
give insights into the SE breakdown and the general effects 
that supercooling has on liquids. 

First, we consider the HS fluid. Figure |6]displays the trans- 
port coefficients and the SE relationship for this system as a 
function of the structural order order parameter, — S2, which 
itself depends only on packing fraction 0c- Although the ba- 
sic scaling relations outlined above describe the behavior of 
attractive fluids along isochores, they apply equally well for 
athermal fluids (such as this) along isotherms, keeping in mind 
that the glass transition in the latter is approached by compres- 
sion and not by cooling. One important aspect to note in Fig.|6] 
is that there is, in fact, a crossover between "equilibrium" (low 
— S2) and "supercooled" (high — S2) states. In fact, Eq. (|6]l 
provides a quantitative fit to the data for packing fractions in 
the range 0.35 < 0c < 0.55 (1.4 < -S2 < 3.9). Interest- 
ingly, the point at which both transport coefficients diverge 
from this equilibrium scaling relationship closely coincides 
with the packing fraction {(f)c ~ 0.55) where D-q — 1.2/(27r), 
i.e., the breakdown of the SE relation. 




FIG. 6: Transport properties as a function of structural order param- 
eter ~S2 for the HS fluid described in the text, (a) Self-diffusivity 
D, (b) viscosity rj, and (c) the SE relationship. In (a) and (b), the 
blue and red lines represent fits to Eq.[6]for the equilibrium and su- 
percooled states, respectively. 



It is also evident from Fig. |6]that one can use the scaling 
form of Eq. ^ with different pairs of coefficients to approx- 
imately describe the S2 dependencies of the transport coeffi- 
cients for the supercooled HS fluid. The fits of Eq. (|6]l for 
these supercooled states, however, are not as accurate as those 
for the equilibrium fluid data below the crossover. In fact, our 
only goal in fitting the supercooled liquid data to this expo- 
nential form is that it allows us to extract simple quantitative 
measures [B'j^ and B[-^] of the couplings that exist between 
the transport coefficients and the static structure of the fluid. 
It can be seen both from the raw data in Fig. |6] and from the 
values of these coefficients for the "equilibrium" and "super- 
cooled" HS fluid [{Bd = 1.3, B,, = -1.4) and (B'jj = 6.1, 
B'^ = —10.1), respectively] that the breakdown of the SE re- 
lation coincides with qualitative change in how static structure 
correlates to D and 77. In the equilibrium fluid, both transport 
coefficients show weak, and roughly equivalent, couplings to 
S2- However, after the SE breakdown, the transport coeffi- 
cients develop much stronger couplings to structural order 
This is presumably due to the integral role that cooperative 
structural rearrangements play in the relaxation of deeply su- 
percooled liquids. 

It is clear from Fig. |6]that the reason that the SE relation 
breaks down in this system is because the viscosity of the 
deeply supercooled HS fluid becomes much more sensitive to 
changes in static structural order than the self-diffusivity. This 
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FIG. 7: (a) Scaled self-diffusivity DT'^^^, (b) scaled viscosity 
and (c) SE relationship D-q/T versus structural order pa- 
rameter ~S2 for the SW-SRA fluid at several packing fractions (jic. 
Self-diffusivity D, viscosity rj, and the SE relationship Drj for the 
HS fluid are provided for comparison. Filled and open symbols rep- 
resent the high (T > 0.5) and low (T < 0.5) temperature branches 
of the SW-SRA fluid, respectively. Dashed lines in (a) and (b) are 
fits of the SW-SRA data to Eq. ^ and the dashed lines in (c) are the 
products of the respective fits in (a) and (b). Arrows in (a) indicate 
the general direction of increasing T. Red and blue lines have the 
same meaning as those in Fig. [6] 



is consistent with the observations of previous studies that 
have correlated the breakdown of the SE relationship to the 
onset of heterogeneous dynamicsj i^i36.37i38,39 Heterogeneous 
dynamics is typically characterized by the presence of many 
particles that transiently exhibit exceedingly high or low val- 
ues mobility relative to the mean. The highly mobile parti- 
cles have been shown to readily diffuse distances on the order 
of a particle diameter by so-called "hopping" motions. It is 
presumably the presence of these highly mobile particles that 
allows the self-diffusivity to maintain a weaker coupling to 
static structure than the viscosity. 

Displayed in Fig. Q are the scaled transport coefficients, 
DT-i/^ and r/T-i/^, of the SW-SRA fluid as a function of 
the structural order parameter — S2. The T~^/^ factor is in- 
cluded here to remove the trivial "thermal velocity" contribu- 
tion to these transport coefficients, which then allows them to 
be directly compared to the dimensionless values of D and 
T] of the HS fluid. Along these lines, it should be noted that 




FIG. 8: Values of the coupling coefficients and — B', from fits 
of the supercooled SW-SRA data to Eq. ^ (shown in Fig.|7j for the 
self-diffusivity D and viscosity rj, respectively. Also shown are the 
(negative) sum — {B'jj + B^) and the ratio ~{B'^/ B'q) of the expo- 
nents. The dashed line represents the HS fluid value for — (B^/B^) 



the high temperature data points for DT and 
Dri/T of the SW-SRA fluid (T > 0.5, denoted in Fig.|2]by 
filled symbols) are indeed approximately described by the cor- 
responding data for D, 77 and D-q of the HS fluid. One conse- 
quence of this is that the breakdown of the SE relation for the 
SW-SRA fluid upon approaching the repulsive glass transition 
by heating occurs at approximately the same value of — S2 as 
the breakdown of the SE relation for the HS fluid upon com- 
pression. In contrast, the value of — S2 of the SW-SRA fluid at 
the breakdown of the SE relation upon cooling toward the at- 
tractive glass is different for each packing fraction studied. As 
might be expected, fluids with lower packing fractions show 
departures from the slip limit of the SE relation for conditions 
where they exhibit lower amounts of translational structural 
order 

Do the SW-SRA transport coefficients also follow "equilib- 
rium" and "supercooled" exponential scaling branches when 
plotted versus — S2? For 0c < 0.55 and intermediate tempera- 
tures, the data approximately collapse onto the same exponen- 
tial scaling relation obtained from the fit of the equilibrium HS 
fluid data (shown in all panels of Fig.|7]as a blue line). At suf- 
ficiently high or low temperatures, however, the SW-SRA data 
transitions to "supercooled" exponential scalings with differ- 
ent coupling coefficients B'jj and B'^. Interestingly, these tran- 
sitions closely coincide with the breakdown of the SE relation. 
In all cases, similar to the HS fluid, D exhibits a considerably 
weaker dependence on — S2 than does rj (i.e., — B,', > B'jj, see 
Fig. [8]l after the SE breakdown. Although the difference be- 
tween the coupling coefficients increases with (pc, the relative 
magnitude remains approximately constant (—B'yS^ ^1.7) 
and very similar to that of the HS fluid. This latter structure- 
property connection is general in the sense that it approxi- 
mately holds for these model liquids as they become "super- 
cooled" via heating, cooling, or compression. 

We would like to conclude this section with a speculation. 
Although it is not yet possible to equilibrate molecular dynam- 
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the HS fluid. Therefore, it is reasonable to suspect that the 
structural order parameter of the high temperature (repulsive) 
glass — S2,i? will be close to that of the HS glass [Fig.|9{a)]. 
At low temperatures, the transport coefficients of the S W-SRA 
fluid as a function of — S2 also appear as if they may asymptote 
to a different, lower, limiting value [Fig.|9ja)], i.e., a value of 
the structural order parameter — S2.A characteristic of the at- 
tractive glass. The resulting picture is that the repulsive and 
attractive glass lines closely follow two "iso-S2" curves, with 
a sharp transition between the two in a narrow temperature 
range, shown schematically in Fig. |9|b)). Since S2 can be 
readily obtained from static pair correlations, this is a specula- 
tion that could be tested via experiments of SRA colloidal flu- 
ids4S It would also be interesting to explore in future studies 
the extent to which liquid-state (e.g., mode-coupling) theories 
for SRA fluids are able to reproduce the empirical connection 
between S2 and dynamics found in our simulations and to test 
the predictions that they make about how the structural order 
varies along the ideal repulsive and attractive glass lines. 



(b) \ 
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FIG. 9: (a) Schematic representing a speculation about how the 
scaled self-diffusivity vs — S2 isochores of the SW-SRA fluid 
[Fig. Ha)] might behave as the repulsive and attractive glass tran- 
sitions are approached. The quantities S2,r and S2,a represent the 
limiting values of S2 for the repulsive and attractive glasses, respec- 
tively, (b) The black and red curve are the proposed iso-S2 loci in the 
T — 4>c plane at the charactertic repulsive and attractive glass values, 
respectively (discussed in the text). The portions of these lines that 
are solid represent the hypothesized glass transition. The yellow line 
is a narrow transition region where the glass line is proposed to cross 
between the iso-S2 curves. 



IV. CONCLUSIONS 

We have presented new molecular simulation data for vis- 
cosity, self-diffusivity, and the local structural ordering of both 
a hard-sphere fluid and a square-well fluid with short-range 
attractions (relative to the particle diameter). We found that 
the latter system has a region of mobility anomalies in the 
temperature-packing fraction plane, where its self-diffusivity 
increases upon isochoric cooling. This region is entirely en- 
closed within a wider set of state points where the fluid's pair 
correlations strengthen upon isochoric heating. This type of 
"cascade of anomalies" is very similar to that found in recent 
simulations of liquid water, and it follows from a broader con- 
nection between static structure and dynamics in condensed 
phase systems. 

Both the hard-sphere and square-well fluids show that the 
breakdown of the Stokes-Einstein relation upon supercooling 
occurs for conditions where viscosity and self-diffusivity de- 
velop different couplings to the degree of pairwise structural 
ordering of the liquid. We discussed how these couplings re- 
flect dynamic heterogeneities. Finally, we provided an ex- 
perimentally testable hypothesis about how repulsive and at- 
tractive glasses may be generally characterized by two dis- 
tinct levels of short-range structural order. In future work, 
we will investigate whether there are similar connections be- 
tween non-equilibrium dynamics (e.g. shear-dependent vis- 
cosity) and structural order in these systems. 
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